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\  Technical  Report  Summary 

..  Ai 

I.  The  great  influence  of  the  way  in  which  inelastic  attenuation  is 


modelled  on  the  results  obtained  for  seismic  source  parameters  from  wave 
spectra  has  led  to  an  investigation  of  better  ways  of  arriving  at  the 
required  corrections  to  raw  seismogram  spectra.  A  procedure,  based  on  a 
relaxation  model  of  attenuation,  for  using  a  frequency-dependent  Q  in  the 
program  for  computing  source  spectra  has  been  developed.  The  resulting 
effective  Q's  tie  together  apparently  disparate  results  obtained  by  others 
for  short-and  long-period  attenuation.  Moments  and  corner  frequencies 
for  12  central  Asian  earthquakes  in  four  selected  seismic  zones  in  western 
China  and  Siberia  have  been  determined  from  P-  and  S-body  wave  spectra.  The 

pi  pc 

range  of  calculated  moments  is  10  to  10  dyne-centimeters.  The  corner 
.  frequencies  from  P-wave  spectra  decrease  with  increasing  moment  in  a  way  that 
is  at  least  qualitatively  as  expected.  However,  the  S-wave  corner  frequencies 
are  found  to  be  independent  of  moment  for  this  small  sample  of  events,  a 
result  that  calls  for  further  consideration  of  the  validity  of  the  analysis. 

2.  The  computation  of  seismic  source  parameters  from  the  properties  of 
the  seismic  signal,  expressed  in  either  the  time-domain  or  frequency  domain, 
depends  on  the  investigator's  choice  of  source  model  from  among  several 
competing  ones.  An  opportunity  to  test  competing  models  empirically  is,  there¬ 
fore,  of  great  value.  The  Haicheng,  China,  earthquake  of  February  4,  1975  is 
the  only  Asian  earthquake  for  which  we  have  detailed  information  about  source 
parameters  from  on-site  geodetic  measurements  and  aftershock  locations  based 
on  a  local  network.  We  are  using  these  local  data  to  calibrate  techniques 
for  determining  the  same  source  properties  from  teleseismic  data.  The  data 
have  been  screened  for  quality  and  the  spectra  are  being  computed.  The 


waveforms  will  be  compared  with  synthetic  seismograms  calculated  on  the  basis 
of  Archambeau's  source  theory. 

3.  Previous  work  under  this  grant  by  Boucher  and  Herrmann  has  led  to 
the  development  of  a  technique  for  estimating  focal  depth  and  focal  mechanism 
of  crustal  events  from  synthesized  long-period  seismograms.  An  attempt  has 
been  made  to  automate  the  Herrmann  technique,  but  the  program  is  not  yet 
operational  because  of  the  existence  of  secondary  maxima  in  the  correlation 
coefficients  used,  such  that  the  program  will  sometimes  find  one  of  these 
and  select  a  wrong  answer. 

The  resolution  and  accuracy  of  the  technique  depends  on  the  adequacy  of 
the  models  selected  for  source  and  receiver  crustal  structures.  A  program 
has  been  developed  that  uses  data  from  earthquakes  recorded  on  a  local 
network  to  simultaneously  locate  the  set  of  earthquakes  and  fix  a  crustal 
model.  The  model  always  incorporates  a  linear  increase  of  velocity  with  depth 
in  the  crust.  The  parameters  calculated  by  the  generalized  inversion  pro¬ 
cedure  are  the  thickness  of  the  crust,  the  surface  velocity,  the.  rate  of  increase 
of  velocity  with  depth  in  the  crust,  and  the  Pn  velocity.  For  those  places 
at  which  the  required  data  are  available,  this  program  will  yield  results  that 
will  improve  significantly  the  output  of  the  focal  parameter  procedure. 


1.  Spectral  Classification  of  Asian  Seismic  Events 
Gary  Lundquist 

Introduction. 

We  have  used  averaged  seismic  source  spectra  to  determine  both  seismic 
source  parameters  and  transmission  properties  of  the  Earth.  Though  each 
estimation  involved  assumptions  about  the  other,  we  feel  that,  by  varying 
parameters  and  drawing  upon  externally  derived  theories,  unknowns  in  both  the 
source  and  transmission  may  be  estimated. 

The  computation  of  the  source  spectra  involved  several  stages.  In  the 
first  stage,  smoothed  autospectra  were  computed  for  the  available  seismograms 
of  12  central  Asian  earthquakes,  3  in  each  of  the  4  geological  provinces 
previously  selected  for  study.  Semi-annual  report  No.  2  of  October,  1975, 
described  the  computation  procedure.  At  that  time,  stations  with  bad  signal- 
to-noise  ratio  were  dropped  from  further  processing,  though  some  borderline 
cases  were  retained.  In  the  second  stage,  initial  corrections  were  made  to  the 
seismogram  spectra,  including  instrumental  response,  geometric  attenuation  (or 
spreading),  crustal  response,  acoustic  impedance  contrast,  and  anelastic 
attenuation.  We  considered  the  correction  for  radiation  pattern  and  concluded 
that  such  a  correction  was  inappropriate  for  crustal  earthquakes.  First,  pP 
arid  sP  were  also  in  the  time  window,  and  we  would  have  been  correcting  the 
energy  transported  by  those  phases  with  a  perhaps  significantly  different  radiation 
pattern  computed  for  direct  P.  In  addition,  we  have  found  accurate  focal 
mechanisms  difficult  to  determine.  The  events  are  small;  station  distribution 
is  not  optimum,  and  the  adopted  take-off-angles  from  the  source  are  almost 
certainly  in  error  due  to  complexities  in  the  crust. 

Once  initial  corrections  had  been  made,  the  next  step  was  to  examine  the 
adequacy  of  some  of  those  corrections.  In  the  case  of  the  spreading  correction, 
we  found  evidence  of  consistent  undercorrection  of  near  stations.  That  is, 


some  travel  paths  in  the  10-25°  range  in  central  Asia  suffer  much  more  rapid 
spreading  than  paths  in  an  average  Earth  model.  The  effect  of  this  under¬ 
correction  is  to  cause  a  spread  in  levels  of  individual  station  spectra,  and 
to  make  the  seismic  moments  slightly  lower  for  those  events  which  rely  on 
near  stations.  The  shape  of  the  averaged  spectra  will  not  be  dominated  by 
one  high  level  station,  however,  since  we  employ  a  logarithmic  average  rather 
than  an  RMS  average.  The  interpretation  of  the  evidence  on  spreading  will 
require  further  study,  and  should  probably  be  combined  with  a  ray  tracing 
experiment. 

The  investigation  of  the  inadequacy  of  standard  anelastic  attenuation  models 
shows  more  immediate  results  in  terms  of  Earth  properties  and  is  the  subject 
of  the  first  major  section  of  this  chapter.  Briefly,  we  found  a  range  of  models 
available,  some  derived  from  surface  waves  and  others  for  high-frequency 
body  waves.  Rather  than  discard  any  model,  we  have  generated  a  modulation  of 
attenuation  based  upon  the  relaxation  theory  of  dissipation  mechanisms  which 
resolves  the  differences  between  models. 

Finally,  we  went  back  to  the  raw  spectra  and  computed  station  spectra 
employing  the  corrections  that  we  felt  were  the  best  available.  Those 
spectra  and  their  interpretation  are  the  subjects  of  the  second  part  of  this 
chapter.  The  consistency  of  magnitude-moment  curves  and  moment-corner  frequency 
graphs  provide  checks  on  the  computation  and  correction  techniques  used. 

The  anelastic  attenuation  of  seismic  body-waves. 

Our  treatment  of  anelastic  attenuation  began  with  the  basic  assumption 
that  the  specific  quality  factor,  Q’1  ,  was  not  frequency- independent. 

This  decision  stemmed  from  the  differences  between  two  California  Institute  of 
Technology  models.  CIT  208  was  derived  from  short-period  body  waves,  and  for 


shallow  earthquakes  could  be  represented  by  a  constant  T/Qg^  ratio  of  ^0.45  , 
where  T  is  travel  time  and  Qg^  is  Q  integrated  over  the  travel  path. 

CIT  11CS2-QM  ,  on  the  other  hand,  was  determined  from  surface  wave  data  and 
gave  T/Qe^^.  n,  1.0  .  The  discrepancy  between  the  T/Qg^  ratios  could  imply 

an  error  in  one  model,  but  we  have  chosen  to  accept  both  as  correct,  and  to  take 
Q"1  is  frequency  dependent. 

We  may  write  the  energy  of  a  harmonic  wave  in  terms  of  its  initial  value, 
Eq(w)  and  a  spreading  factor  G(x)  ,  as  (Anderson,  1967) 
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where  k  is  a  complex  wave  number,  <  =  <'  +  i<"  .  The  specific  quality 

factor,  Q"1  ,  is  defined  as  the  fraction  of  energy  dissipated  per  wavelength. 
If  the  energy  loss  is  small,  then  to  first  order  the  energy  carried  is  given 
by  (1.1).  Neglecting  common  factors, 
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For  small  ’  ,  the  series  expansion  of  an  exponential  gives 
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We  may  distinguish  between  P-wave  and  S-wave  attenuation  by  writing  the 


respective  wave  numbers  in  terms  of  complex  velocities. 
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where  k  and  p  are  complex  bulk  and  shear  moduli,  respectively.  For  the 
case  of  small  attenuation, 
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(1.5) 
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Equations  (1.5)  give  the  explicit  dependence  of  anelastic  attenuation  on 
the  imaginary  parts  of  the  elastic  moduli.  It  is  important  to  note  that  both 
P  and  S  waves  suffer  attenuation  from  shear  mechanisms  (u"  ?  0),  but  P 


waves  may  suffer  additional  attenuation  from  bulk  losses  (k’‘  t  0).  In 
■1  , 


assuming  Q  is  a  function  of  frequency,  we  are  assuming  that  p11  and  k1 
are  frequency  dependent.  Bulk  losses  are  generally  assumed  to  be  small,  and 
setting  k'1  =  0  in  (1.5)  gives 
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so  that. 


=  2.25  for  Poisson's  ratio  =  0.25.  If 


(1.6)  implies  q 
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=1.0  ,  then 


'u  4.0 


We  will  relate  these  definitions  of  Q  to  relaxation  phenomena  through 
the  "standard  linear  elastic  solid"  model  (Mason,  1958).  The  model  may  be 
visualized  as  a  spring  of  stiffness,  Mj  ,  in  series  with  a  parallel 
combination  of  dashpot  (viscosity  02  )  and  spring  (stiffness  M2  ).  At  high 
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frequencies,  or  upon  sudden  application  of  a  stress,  the  system  has  an 
instantaneous  elastic  response  controlled  by  Mj  ,  the  unrelaxed  system 
elastic  modulus.  At  very  low  frequencies,  or  upon  application  of  a  constant 
stress,  the  system  modulus  is  the  relaxed  stiffness  of  the  series  springs, 

M  -"I"? 

MR  ~  Mj+M2 


The  anelastic  dashpot  controls  the  frequency  range  over  which  the  change  in 
modulus  goes  from  to  MR  ,  and  therefore  the  frequency  band  over  which 
absorption  takes  place. 

For  a  harmonic  stress,  the  stress-strain  relation  is  (Liu  et  a_l_. ,  1976;. 

a  +  Ta  6  =  MR(e*T£  t)  (1.7) 


where  xa  =  (M^  +  M2)  is  the  stress  relaxation  time  for  constant  strain, 

and  t£  =  n2/M2  is  the  strain  relaxation  time  for  constant  stress.  Since 
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Thus  we  obtain  the  complex  modulus  as 
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where  t  =  /r’a ’t£  is  the  average  system  relaxation  time,  and  C^,  Mg)  is 
a  constant  depending  upon  the  elastic  moduli  of  the  system.  In  (1.10)  we 
have  isolated  the  frequency  dependent  behavior  inside  i.ie  square  brackets. 

That  factor,  which  we  call  R(u;)  ,  gives  a  bell -shaped  graph  of  Q'1  vs. 

frequency. 

The  above  are  all  well  known  results,  and  any  medium  which  may  be  modelled 
as  a  standard  linear  solid  will  show  absorption  over  a  limited  frequency  band 
and  appear  elastic  at  other  frequencies,  though  different  elastic  moduli 
apply.  The  Earth  is  much  more  complex  than  this  standard  solid  however, 
and  such  a  model  applies  only  over  limited  frequency  ranges.  Liu  _et  aj[,  (1976), 
however,  believe  the  model  can  be  applied  about  relaxation  times  so  closely 
spaced  that  the  various  absorption  bands  overlap,  and  Q"^  is  essentially 
constant.  This  requires  frequent  changes  in  the  imaginary  elastic  moduli, 
and  many  mechanisms  have  been  described  which  will  give  the  required  spectrum 
of  relaxation  time,  and  an  essentially  frequency  independent  Q~* 
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The  application  of  anelastic  attenuation  correction  based  on  the  constant 

T/Q  ^  ratios  implied  by  (1.6),  however,  leads  to  inconsistencies. 

Conservation  of  energy  requires  that  amplitude  spectra  must  decay  at  least  as 
-3/2 

fast  as  u  at  high  frequencies.  And  interference  under  Huyghen's 
principle  indicates  that  the  decay  should  occur  when  wavelengths  are  on  the 
order  of  fault  size  or  smaller,  implying  a  corner  frequency  well  within  the 
passband  of  WWSSN  seismographs.  Also,  the  shortest  time  constant  of  rupture 
is  probably  on  the  order  of  the  rise  time  or  smaller,  implying  smooth  or 
continuous  acceleration  to  frequencies  greater  than  the  inverse  of  the  rise 
time.  Those  frequencies  should  also  be  in  the  passband  of  WWSSN  short-perioa 

_3 

instruments,  and  should  decay  at  least  as  fast  as  u>  .  Observed  slopes 
at  high  frequencies  after  initial  corrections  were  positive,  rather  than 
negative,  especially  for  S-wave  spectra,  implying  an  error  in  the  correction 
technique. 

The  total  anelastic  attenuation  correction  is  give-i  by 


uiT 

2Qeff 


(1.11) 


The  slope  of  the  corrected  spectrum  is  affected  by  either  changing  T/Q  ^ 
or  making  a  function  of  frequency.  We  have  chosen  the  latter,  and 
modulate  Q~*  according  to 


Q_1(z,  u)  =  Q_1(z)  ,  *  1 

Q_1  (z ,  w)  =  Q_1(z)  R(c)  ,  w  >  4- 
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where  R(u>)  is  defined  in  (1.10).  One  aspect  of  the  substitution  of  Rl.  ' 

into  the  exponent  of  (1.11)  is  the  fact  that  is  an  integration  of  Q  »  « 

over  the  travel  path.  According  to  (1.12),  we  have  assumed  R(„)  to  be 
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of  path,  so  that  it  comes  out  of  the  integration.  This  is  obviously  a  simplifi¬ 
cation,  and  the  property  of  path  independence  must  be  tested.  In  particular, 
we  expect  a  different  modulation  for  deep  events. 

The  result  of  the  modulation  defined  by  (1.12)  is  that  we  model  a  hole  in 
the  spectrum  of  relaxation  times.  The  decay  in  Q"*  has  the  shape  of  a 
similar  decay  in  a  linear  elastic  solid,  with  the  beginning  of  the  decay  at 
w  =  1/t  .  Figure  1.1  shows  the  total  anelastic  attenuation  correction  for 

P  and  S  waves  as  a  function  of  x  ,  where  the  unmodulated  correction 
assumes  constant  T/Qeff  ratios  of  1.0  and  4.0  for  P  and  S  waves,  respectively. 
In  Figures  1.2  to  1.5  we  show  spectra  corrected  with  various  x  modulations 
for  each  of  our  four  geographic  areas. 

The  P-wave  spectra  yield  to  simple  interpretations,  x  =  0.3  is  adequate 

o. 

for  all  12  events  studied  thus  far,  though  events  dominated  by  near  stations 
could  just  as  well  use  larger  x&  .  The  choice  of  x^  is  made  according 
to  the  slope  criterion  previously  discussed,  and  an  u  decay  rate  is  sought. 
The  choice  of  xQ  =  0.3  gives  an  acting  T/Qef.^  ratio  of  1.0  at  20  sec 
period  and  0.5  at  .5  sec  period,  effectively  combining  the  two  CIT  models 
using  a  standard  elastic  solid  decay  rate,  x  is  also  path  independent,  so 
far  as  this  study  is  concerned. 

The  interpretation  of  the  S-wave  spectra  is  not  so  simple.  In  this  report 
we  present  the  results  of  using  T/Q  ^  =  4.0  ,  and  comment  on  the  conclusions 
which  would  have  to  be  made  if  this  model  were  correct.  If  the  choice  of 
Xg  were  made  on  the  basis  of  high  frequency  decay,  x^  =  0.7  would  be  the 
minimum  value.  However,  a  subjective  estimate  of  spectral  quality  leads 

this  investigator  to  prefer  xQ  =  1.1  .  For  x„  <  1.1  ,  most  of  the  S- 

p  p 

wave  spectra  had  definitely  peaked  character,  and  the  interpretation  of  a 
long  period  level  for  seismic  moment  is  difficult. 
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In  either  case,  r0  <  0.7  does  not  give  an  acceptable  slope,  implying 

P 

that  P  waves  are  attenuated  to  higher  frequencies  than  are  S  waves.  That 
is,  the  spectrum  of  relaxations  which  are  controlled  by  the  shear  modulus 
has  a  decay  toward  high  frequencies  from  a  peak  at  about  .15  Hz,  while  a 
separate  dissipation  mechanism  operating  on  the  bulk  modulus  maintains 
p-wave  attenuation  out  to  .5  Hz.  The  half  amplitude  points  for  Q"1  and 
Q"1  are  at  about  0.5  and  2  Hz,  respectively,  for  Tg  =  1.1  and  =  0.3 
sec.  If  t0  =  0.7  ,  the  peak  would  be  at  .2  Hz  and  the  half  amplitude  at 

P 

.7  Hz.  The  range  of  acceptable  ig  never  permits  a  complete  overlap  of 
the  P  and  S  attenuation  spectra.  Figure  1.6  diagrams  these  results  on  a 
log- log  graph. 

Two  bulk  loss  mechanisms  have  been  studied  in  the  literature.  Biot 
(1956,  1956a)  studied  the  propagation  of  stress  waves  in  a  saturated 
porous  solid.  In  his  model,  the  pores  were  interconnected,  and  a  compressible 
fluid  was  allowed  to  move  with  respect  to  the  solid  matrix.  His  results  show 
peaked  attenuation,  with  Q  decaying  as  u  2  at  high  frequencies.  The 
location  of  the  peak  depends  in  a  complex  way  upon  fluid  viscosity,  velocity 
and  density,  the  frictional  forces  on  the  fluid,  the  pore  (tube)  diameter 
and  length  and  the  mathematical  model  of  flow.  There  is  sufficient  flexibility 
in  his  parameterizations  that  such  a  mechanism  could  be  called  upon  to  explain 
the  observed  P  wave  attenuation  peak.  However,  the  decay  as  only  u. 
would  require  a  different  choice  of  t  ,  a  test  which  has  not  been  made. 

The  second  possible  bulk  loss  mechanism  was  examined  by  V'aisnys  (1968). 

He  proposes  that  a  material  which  is  in  equilibrium  very  near  its  melting 
point  may  change  phase  under  stress.  And  if  the  stress  is  harmonic,  the  entropy 
production  would  vary  with  frequency,  being  peaked  for  frequencies  about  1/-.  , 
che  characteristic  rate  constant  of  the  process.  Essentially,  the  medium  would 
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FI  pure  1.6:  The  frepuency  dependent  corrections  tc  the  r.ei  sr:o''r\  spectra.  At  th" 
hotter  are  the  instrument  responses  for  IWSSM  seisi-orranhs.  At  the  to" 
n-1  is  olotteri  vs.  freoucncy  for  ■  ’ .  1  a  no  0.3. 


tend  to  freeze  during  one  half  cycle,  then  thaw  during  the  other  half  cycle. 
Since  more  energy  is  required  to  order  a  lattice  than  to  disorder  a  lattice, 
the  effect  would  be  a  slightly  greater  fraction  of  melt  after  the  wave  has 
passed.  Since  this  mechanism  requires  pressure  changes,  it  would  attenuate 
only  P  waves. 

Vaisnys  examines  three  cases:  (1)  an  ideal  solution  partial  melt; 

(2)  solid  in  equilibrium  with  its  own  melt;  and  (3)  solid-solid  phase 
transformation.  For  each  he  finds  attenuation  to  be  peaked  in  exactly  the 
form  we  have  described  for  a  standard  linear  solid.  The  relaxation  times  are 
10  sec,  90  sec  and  greater  than  10^  sec,  respecti vely.  Thus  all  time  constants 
are  greater  than  =  0.3  as  observed  here-  However,  though  the  thermodynamic 
character  of  the  medium  is  completely  parameterized,  Vaisnys  makes  no  effort 
to  compare  those  properties  to  upper  mantle  properties.  Thus  this  mechanism 
might  also  apply. 

We  state  again,  that  the  separation  of  distinct  P  and  S  attenuation 
mechanisms  relies  upon  both  our  choice  of  t  and  T/Qeff  .  We  will  see  in 
the  next  section  that  our  choice  of  parameters  for  S  is  likely  in  error,  since 
corner  frequencies  tend  to  concentrate  around  .25-. 3  Hz,  independent  of  the 
size  of  the  events.  Also,  while  P  wave  spectra  sometimes  appear  to  be  two 
cornered,  the  S  wave  spectra  rarely  are.  Thus  there  is  good  reason  to  question 
the  validity  of  the  choice  of  parameters  for  S  and,  perhaps,  of  the  whole 
concept.  This  writer  feels  the  beauty  of  the  relaxation  concept  of  specific 
dissipation  mechanisms  as  well  as  the  consistency  of  the  P  wave  results 
justify  the  present  effort,  as  well  as  further  work. 

Presentation  of  spectra. 

Average  seismic  body-wave  spectra  are  presented  in  Figures  1.7  to  1.18, 
with  three  earthquakes  studies  in  each  of  (1)  the  Tien  Shan  fold  belt;  (2)  the 
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rigure  1.7:  Average  source  snectrur,  Tien  Shan  event  of  72-04-09.  SPP  from  7, 
LPP  from  7. 
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.9:  Average  source  s^ectrun,  Tien  Shan  event  of  69-0?- 1 1 .  SPP  fro 
LPP  frnn  7,  SPt.  froro  3,  IPS  fron  T>. 
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Finure  1.16:  Average  source  s^ectrur',  l.-ikc  H  i ka  1  event  o'  7P-0*'-09. 
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;e tan  plateau,  (3)  the  block-faulted  Tsaidam  basin  and  (4)  the  Lake  Baikal 
rift  zone.  Separate  and  independent  estimates  were  made  of  short-period 
P  and  SH  spectra  and  long-period  P  and  SH  spectra.  In  each  case,  the  best 
available  data  were  used,  whether  or  not  corresponding  data  was  available  for 
the  other  components. 

A  crude  estimate  of  the  reliability  of  the  spectra  is  given  by  the  number 
of  stations  in  each  average,  as  noted  in  the  figure  caption.  We  consider 
an  average  of  six  or  more  stations  to  be  good,  while  an  average  of  four  or 
less  is  questionable.  Neither  Tsaidam  events  of  72-08-30  has  sufficient  short- 
period  S  (SPS)  data,  and  those  spectra  are  presented  only  for  comparison. 

As  mentioned  in  the  introduction,  near  stations  are  systematically 
undercorrected  for  geometric  attenuation.  That  is,  averages  dominated  by 
near  stations  tend  to  be  low.  The  Tsaidam  events  of  72-08-30  should  have 
i'.ightly  low  seismic  moments,  and,  in  particular,  the  SPS  spectra  are  low 
by  a  factor  of  about  2  to  5.  The  corner  frequencies  were  picked  on  the  LPS 
spectra. 

The  spectra  are  plotted  to  seismic  moment  scale,  and  the  choice  of 
moment  is  shown  on  each  figure  by  an  arrow,  and  listed  in  Table  1.1.  Moments 
are  plotted  vs.  earthquake  magnitude  in  Figure  1.19,  where  the  horizontal 
bars  show  the  difference  between  body-wave  magnitude,  m^,  and  surface  wave 
magnitude,  M  .  The  P-wave  data  fit 

M  =  -1.7  +  1.3  log  MQ 

Kanamori  and  Anderson  (1975)  find  a  slope  of  2/3,  which  we  cannot  possibly 
fit.  The  data  presented  here  would,  however,  fall  into  the  scatter  of 
'  amori  and  Anderson's  points  at  the  low  magnitude  end. 
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TABLE  1.1 


I - f4 


Macmitude 


Seismic  Moment  Corner  Frequency 


1 

mb 

Ms 

\ 

Mo 

s 

f% 

(Hz) 

Tien  Shan 

, 

A 

72-04-09 

5.73 

5.43 

4.0x1a24 

.90 

■1 

70-06-05 

6.17 

6.67 

l.OxlO26 

3.2xl025 

.26 

.26 

i 

s 

69-02-11 

6.15 

6.38 

3.6xl025 

2. OxlO25 

.43 

.27 

- 

Tibet 

•  ; 

71-05-22 

5.85 

6.04 

25 

1.0x10“ 

9. lxlO24 

.42 

.28 

r 

? 

71-04-03 

5.91 

5.99 

6.8xl024 

1.2xl025 

.59 

.27 

.  4 

67-08-15 

6.17 

6.67 

1.9xl024 

1.8xl024 

.66 

.31 

i, 

Tsai  dam 

r 

72-08-30:18 

5.22 

4.99 

1.5xl024 

4.0x1 O23 

.60 

.30 

72-08-30:15 

5.26 

4.86 

2AxlOU 

9. lxlO23 

.51 

.30 

\ 

63-04-19 

6.10 

6.36 

7.2xl025 

1.4xl026 

.37 

.16 

V  ' 

s 

Lake  Baikal 

72-08-09 

4.7 

5. OxlO23 

1.2 

70-08-13 

4.7 

2.8xl023 

.86 

» 

■  V 

p 

67-02-11 

5.15 

5.16 

9. lxlO23 

.80 

*  v 
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Note  that  seismic  moments  computed  from  S-wave  spectra  are  smaller  than 
those  computed  from  P-wave  spectra  with  the  major  exceptions  of  the  Tsaidam 
event  of  63-04-19  and  the  Tibet  event  of  71-04-03.  The  levels  are  most 
consistent  in  Tibet. 

Corner  frequencies,  f  ,  are  also  noted  by  arrows  on  spectral  plots, 
listed  in  Table  1.1  and  plotted  in  Figure  1.20  on  an  M  -fQ  diagram.  fQ 
is  defined  as  the  intersection  of  a  horizontal  line  through  the  low-frequency 
trend  with  a  line  through  the  high-frequency  asymptote.  Rather  than  assume 

_3 

a  particular  high-frequency  decay  rate,  such  as  ,  and  fit  that  slope 

to  the  data  as  well  as  possible,  a  best  fit  line  was  drawn  through  the  data. 
And  rather  than  average  through  the  entire  short-period  spectrum,  the  quality 
was  judged  with  respect  to  the  passband  of  the  instrument.  The  instrument 
responses  are  presented  in  Figure  1.6  to  demonstrate  that  the  high-frequency 
slope  should  be  estimated  from  frequencies  above  0.5  Hz. 

-3 

If  two  choices  of  decay  rate  could  be  made,  the  slope  closer  to  w 
was  chosen.  The  problem  is  to  distinguish  between  true  slope  changes  and 
holes  in  the  spectra.  Such  questions  arose  on  the  P-spectra  of  71-05-22, 
71-04-03,  and  69-02-11. 

An  examination  of  the  spectra  will  show  that  more  than  half  of  the  P 
spectra  could  be  interpreted  as  having  more  than  one  corner  frequency,  while 
seven  of  eight  S  spectra  have  virtually  the  same  shape.  And  while  P  corner 
frequencies  vary  with  moment,  the  S  corner  frequencies  are  essentially 
constant.  This  inconsistency  between  P  and  S  could  possibly  be  argued  in 
terms  of  source  properties,  P-coda  contamination  of  the  S  wave,  PP  and  PP;' 
contamination  of  the  P  wave  or  other  factors,  but  it  more  likely  represents 
an  insufficient  knowledge  of  the  anelastic  attenuation  operator. 
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Twelve  earthquakes  is  hardly  a  statistically  significant  sample  of  central 
Asian  seismicity,  and  these  first  efforts  do  not  lead  to  a  consistent 
scaling  law  relationship.  If  anything,  these  events  demonstrate  the  expected 
heterogeneity  of  intraplate  earthquakes.  More  data  are  required  to  resolve- 
regional  characteristics. 
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2 .  Focal  Parameters  from  Long-Period  Body  Waves 
W.  Gawthrop 

Synthesis  of  Seismograms. 

Body  wave  synthesis  has  been  used  in  determining  focal  parameters  of  a 
few  crustal  earthquakes  in  Asia.  Focal  depths  and  focal  mechanisms  of  earth¬ 
quakes  also  being  studied  by  analysis  of  their  spectra  (Lundquist,  this  issue) 
have  been  the  object  of  this  study.  Several  practical  problems  have,  however, 
slowed  progress  in  determining  these  focal  parameters.  The  technique  of 
digitizing  long  period  P-waves  in  the  presence  of  high  noise  levels  from 
WWSSN  stations  to  obtain  both  usable  spectra  as  well  as  input  to  body  wave 
synthesis  programs  has  proved  to  be  less  than  successful.  Also,  the  new 
program  to  automatically  determine  focal  parameters  by  comparing  synthetics 
digitized  seismograms  has  not  worked  properly  to  date.  The  new  program 
searches  several  trial  depths  and  focai  mechanisms,  then  iterates  from  the 
best  trial  point  to  obtain  the  maximum  correlation  coefficient.  Local 
maxima  of  correlation,  which  exist  in  the  scattering  of  trial  focal  para¬ 
meters,  confuse  the  program.  Until  some  control  is  put  into  the  new  program, 
data  must  be  of  exceptionally  high  quality,  with  a  high  signal -to-noise 
ratio,  for  a  good  solution  to  be  obtained.  In  addition,  a  means  of  automatical 
estimating  the  errors  in  the  focal  parameters  still  needs  to  be  developed. 

The  focal  parameters  obtained  from  studying  seismograms  in  the  time 
domain  through  careful  analysis  are  in  decreasing  order  of  resol vabi 1 i ty , 
the  depth,  orientation  of  the  fault  planes,  seismic  moment,  duration  of 
faulting,  and  direction  of  propagation  of  rupture.  Figure  2.1  shows  the 
first  20  seconds  of  seven  well-distributed  seismograms  from  a  thrust  ea^thquaxe 
in  the  Tien  Shan  province  of  northwestern  China  on  June  5,  1970.  Relative 


arrival-times  of  the  P,  pP,  ana  sP  arrivals  determine  a  source  depth  of 


MAT 


IS  kilometers,  assuming  the  earthquake  to  be  a  point  source  and  the  crust  to  be 
of  uniform  velocity,  6.0  km/sec.  Since  the  shapes  of  all  seven  seismograms  are 
quite  similar  for  this  event,  the  focal  mechanism  can  only  be  determined  by  care¬ 
ful  determination  of  the  relative  amplitudes  of  the  peaks  corresponding  sc  direct 

and  reflected  wave  arrivals  of  maximum  energy  from  the  source  pulse.  The  area 

26 

under  the  source  pulse  gives  a  seismic  moment  of  2  x  10  dyne  cm. 

The  duration  of  the  pulse  is  about  5  seconds,  with  stations  to  the  east 
requiring  a  slightly  shorter  pulse.  This  shortening  of  the  pulse  to  the  east 
is  probably  due  to  a  doppler  shift  from  an  eastward  rupturing  fault. 

Only  about  15  seconds  of  the  seismograms  can  now  be  accurately  matches. 

This  is  because  most  of  the  early  energy  comes  from  the  P,  pP,  and  sP  waves 
which  have  very  little  dependence  on  the  crustal  structures  at  the  source 
and  receiver.  Synthesis  of  wave  shapes  beyond  about  15  seconds  can  only 
be  considered  reliable  if  we  know  with  confidence,  by  independent  means, 
the  crustal  structures  in  the  areas  of  interest. 

Crustal  Structure  Determination. 

A  program  has  been  developed  to  independently  determine  a  crustal 
velocity  structure  and  Pn  velocity  in  seismically  active  regions  where  locd 
networks  of  seismographs  exist.  The  program  utilizes  a  generalized  inversim 
procedure  to  simultaneously  determine,  from  P-wave  arrival  times,  hypocente- 
locations  of  up  to  12  earthquakes  and  the  velocity  structure  in  the  local 
region.  A  crustal  structure  with  velocity  increasing  linearly  with  deptn  dov.: 
to  a  constant  velocity  upper  mantle  permits  four  accurately  determined 
parameters  to  completely  describe  the  structure.  These  parameters  are 
velocity  at  the  surface,  the  velocity  gradient  with  depth,  the  deptn  of  t-c 
crust,  and  the  velocity  at  the  top  of  the  mantle. 


Individual  data  are  weighted  according  to  tne  observational  error  to 
give  the  data  matrix  unit  variance  for  the  error  analysis.  Differently 
dimensioned  parameters  are  scaled  to  give  each  one  equal  weight  in  the  least 
squares  procedure.  A  Lagrangian  multiplier  is  incorporated  to  increase 
stability  during  successive  iterations.  The  solution  to  the  equations  are 
given  in  matrix  form  as 

x  *  (S(ST  AT  WT  W  A  S  +  e2I)_1  ST  A'  WT  wJb  =  GB 

RESOLUTION  MATRIX  =  GA 

COVARIANCE  MATRIX  =  G(varB)GT  =  GGT 

A 

where  x  is  the  calculated  change  in  the  model  parameters  during  each  iteration, 
S  is  the  diagonal  scaling  matrix,  A  is  the  matrix  of  derivatives  with 
respect  to  the  model  parameters,  W  is  the  diagonal  weighting  matrix,  s 
is  the  Lagrangian  multiplier,  B  is  the  data  matrix,  and  G  is  the  generalized 
inversion  matrix.  Iteration  continues  until  a  convergence  criterion  is 
satisfied. 

Preliminary  tests  in  California  and  Missouri  indicate  that  the  deptn 
of  the  crust  can  be  determine  to  about  25%  and  the  velocity  structure  of  the 
crust  to  about  10%.  The  accuracy  improves  significantly  when  known  sources, 
such  as  explosions,  are  used,  since  in  many  cases,  a  tradeoff  exists  between 
the  focal  depths  of  the  seismic  event  and  the  thickness  of  tne  crust.  It 
is  anticipated  that  the  addition  of  S-wave  arrival  times  to  the  data  set 
will  serve  to  fix  the  focal  depths  and  improve  the  resulting  determinations  of 
crustal  thickness.  It  is  hoped  that  the  application  of  this  program  will 
improve  our  understanding  of  crustal  structure  so  we  can  better  model  long- 
period  seismograms  for  determining  focal  parameters. 


3.  The  Haicheng  Earthquake 


Max  Wyss 

The  Haicheng  earthquake,  =  7.4,  occurred  on  February  4,  1975.  t 
was  located  in  northeastern  China  (latitude  40. 6N,  longitude  122. 6E). 
earthquake  is  at  the  present  time  the  only  Asian  event  for  which  we  have 
detailed  information  on  source  parameters  derived  from  geodetic  and  after¬ 
shock  studies  in  the  source  area.  For  this  reason  this  event  is  ideal  for 
testing  and  calibrating  our  methods  for  deriving  source  parameters  for  Asian 
earthquakes.  To  perform  this  test  we  will  match  the  time  series  and  soect**a 
of  P-waves,  S-waves  and  surface  waves  with  time  series  and  spectra  of  theoretical 
source  models. 

The  Haicheng  earthquake  is  the  famous  event  predicted  oy  Chinese 
seismologists.  Thanks  to  the  prediction  effort,  geodetic  measurements 
were  carried  out  in  the  source  region  before  the  earthquake.  From  the  repetition 
of  these  measurements  and  aftershocks  studies  the  Chinese  derived  the  source 
parameters  as  (Muller,  1976): 

Length  =  70  km 

Width  =  17  km  (0  to  17  km  depth) 

Dislocation  =  55  cm 
Stress  drop  -  5  bars 

The  sense  of  faulting  was  left  lateral  \long  a  plane  striking  N70W,  with  the 
propagation  mostly  towards  the  northwest. 

The  first  step  was  to  screen  the  WWSSN  data  for  seismograms  with  ciea> 
signals  at  distances  suitable  for  analysis.  The  second  step  was  to  digitize 
these  signals.  These  two  tasks  are  now  almost  completed.  Table  3.1  lists  the 
station  for  which  seismograms  have  been  digitized  to  date.  The  first  column 
gives  the  station  abbreviation  (WWSSN).  The  columns  headed  P,  S,  9  and  G  give 
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the  number  of  components  digitized  at  each  station  for  the  respective  phases. 
The  number  2  in  the  P  column  means  that  a  long  period  and  a  short  period 
record  (usually  both  z)  have  been  digitized  at  that  station  for  the  P  analysis. 
A  number  4  in  the  S  column  means  that  all  four  horizontal  components  (long 
and  short  period)  have  been  digitized  for  the  S-phase.  For  the  surface 
waves  (R  and  G)  only  long  period  records  were  digitized. 

At  the  present  we  are  in  the  process  of  finishing  the  third  step  in 
the  analysis.  The  Fourier  transformation  of  appropriate  windows  for  all  the 
data  digitized.  The  fourth  and  final  step  is  just  being  started.  We  are 
generating  theoretical  spectra  and  time  series  with  the  Archambeau  model 
using  the  above  source  parameters  as  input.  The  main  variable  will  be 
the  rupture  velocity.  The  objective  is  to  come  up  with  a  source  model  that 
agrees  with  the  field  observations  and  that  generates  the  *ar  field  seismic 
waves  observed.  We  expect  that  the  results  of  this  study  will  be  ready  for 
publication  by  the  end  of  the  year. 
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TABLE  3.1 


Ha i cheng  earthquake 


Station 

ADE 

ANP 

ARE 

AQU 

BAG 

BHP 

B02 

BUL 

CHG 

COL 

COR 

CTA 

OAL 

ESK 

GIE 

GUA 

HNR 

1ST 

JCT 

KBS 

KEV 

KIP 

LEM 

MSH 

NIL 

NUR 

SHA 

SNG 

STU 
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TUC 

VAL 
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stations 
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components 
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E 
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R 

o 

CD 

G 
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deg 
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1 

77 
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1 

15 
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1 
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31 
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1 
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28 
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O 

c 

82 
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74 

331 

40 

1 

2 

130 

45 

327 
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61 
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1 
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67 
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56 

1 

99 

36 
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1 

54 
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1 

2 

69 
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49 
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